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Abstract 


The attenuation of the intra-band intensity of a snperdeformed band which 
results from mixing with normally deformed confignrations is calculated using 
reaction theory. It is found that the sharp increase of the attenuation is mostly 
dne to the tnnnelling through a spin dependent barrier and not to the chaotic 
natnre of the normally deformed states. 

It is now well established that the intensities of E2 gamma transitions within a snperde¬ 
formed (SD) rotational band show cascades down to low angnlar momentnm M- These 
cascades exhibit the distinct featnre that the intensity remains constant nntil a certain spin 
is reached where-after the intensity drops to zero within a few transitions. The sharp drop in 
intensity is commonly referred to as the decay ont of a snperdeformed band and is believed 
to arise from mixing of the SD states with normally deformed (ND) states of identical spin. 

The earliest theoretical work to implement snch an interpretation ^-IT 
model of the conpling between the SD and ND states. More recently, Refs. 


nsed a statistical 
T^JT^ nsed 


a framework originally developed for the stndy of componnd nnclear reactions to derive 
formnlae for the intensity in a more rigorons fashion (the expressions for the intensity in 
Refs. |l8|-[TI| are dednced from probability argnments). Ref. [^] conclnded that Refs. M 
are valid in the non-overlapping resonance region. Refs. HH fnrther calcnlate the spin 
dependence of the relevant parameters (the electromagnetic widths of the SD and ND states, 
the level density of the ND states and spin dependence of the barrier separating the SD and 
ND wells) which Refs. do not. Two featnres common to Refs. are (i) the nse 

of the Ganssian Orthogonal Ensemble (GOE) to simulate the ND states (ii) the use of the 
“golden rule” to extract a width for the the SD states due to mixing with the ND states. 

Here, as in Refs. [0,|1^] we exploit the similarity between the decay out of superdeformed 
bands and compound nuclear reactions to write the intensity as the sum of average and 
fluctuation contributions. However we use an energy average in place of the ensemble average 
used in Refs. The energy average approach allows the inclusion of the following 


features which are more difficult to incorporate into an ensemble average. 

(i) A hierarchy of complexity in the ND spectrum may be introduced. 

(ii) A statistical model different from the GOE may be used to simulate the ND states, as 
was proposed in Refs. 


14.15 


(hi) A width for the SD states due to mixing with the ND states arises naturally without 
appealing to the “golden rule” whose range of validity has been found to be restricted ||l6|,pT . 
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In Figure ^ we show a schematic plot of the energy of ND and SD bands as a function 
of spin. The observable in which attenuation is seen is the total intensity of two consecutive 
E2 photons in the cascade down the SD band. Let | J) denote an SD conhguration with spin 
J. The relative intensity of the two step transition | J + 2) ^ | J) 3- | J — 2) (relative to the 
intensity of the same two step transition in the absence of mixing with other conhgurations) 
is given by 

1 POO POO 

= / dE,, / dE,,\{J -2\T{Ej_^ + E,,)\J + 2)\H{E - Ej_^- E,,) 

ZUl j_|_2 —oo j—oo 

1 

= — dE\{J-2\T{E)\J + 2)\^, (1) 

where E = F^j +2 — E^^ = Ej takes account of the Hamiltonian of the electromagnetic held, 
Ej being the energy of \J) and E^^ and E^^ the energies of the two consecutive photons. 
The electromagnetic width of | J + 2) is rj _,_2 making 27rrj_,_2 the intensity when there is no 
mixing with the ND states and thus no hux loss from the SD band. Note that in Eq. (|^ 
we ignore the widths of the initial and hnal states for the purpose of calculating the relative 
intensity. 

The transition amplitude is given by 

(J - 2|T(E)| J + 2) = 7j+ 2(J|G(E)| J)7j. (2) 

Here 7 j +2 is the electromagnetic decay amplitude of | J + 2) dehned such that rj _,_2 = 7 j+ 2 ; 
7 j and T} = 7 } are the corresponding decay amplitude and width of \J) whilst the Green’s 
function is given by 

G{E) = {E-H)-\ (3) 

The total nuclear Hamiltonian, which takes the coupling to the electromagnetic held into 
account, is denoted by H. 

The projected Green’s function {J\G{E)\J) may be expressed in terms of its Lorentzian 
energy average {J\G^'''{E)\J) = {J\G{E + y)]J) (energy averaging interval I) plus a huc- 
tuation part 

(J|G|J) = (JlCjJ) + (4) 

where by dehnition the energy average of | J) is zero. Thus Eq. (P for the relative 

intensity may be written 

Fj = F7- + Ffl--, (5) 


where 


and 


FT 

FT = ^ j_JE\(J\G"V)\^ 


( 6 ) 

( 7 ) 
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In this paper we focus our discussion on Ff''. It can be shown that 


{J\G^-{E)\J) 


1 

E-Ej + iTy2-Wjj{E)' 


(8) 


The derivation of Eq. (||) for J) and an expression for | J) using projection 

operator techniques will be reported in a subsequent paper. 

The form of Wjj{E) depends on the specihc model for the Hamiltonian which is em¬ 
ployed. It is our aim to study whether or not the chaotic nature (as classified by random 
matrix theory [RMT]) of the ND states is decisive in explaining the observed attenuation. 
In order to isolate the statistical aspects of the calculation we use two different models dis¬ 
tinguished by whether the tunnelling interaction mixes \J) randomly with the ND states 
{model A) or whether it couples more strongly to certain ND states than others {model B). 
In the latter we shall make the most extreme assumption that \J) couples to only one ND 
state. 

Model A is represented by the matrix 


f Ej Vjn 

y Tjn E^dfi'n 


2 1 , 0 Tj,5n'n J ’ 


n = 1, 


(9) 


where E^ denotes the energies of the N ND states with which \J) mixes due to the real 
tunnelling interaction Vj„. Here T^ is an electromagnetic width which we assume to be 
common to the ND states. With these dehnitions Wjj{E) becomes 


N 


wjj(E )= 


K 


Jn\ 


n=l 


E-E„ + i(rl, + I)/2 


( 10 ) 


The energies En are constructed using the deformed Gaussian orthogonal ensemble 
(DGOE) The DGOE allows a smooth interpolation from Poisson to GOE statistics by 
varying a mixing parameter A from 0 to 1. Thus the En are the eigenvalues of a random 
Hamiltonian h which is real symmetric and whose matrix elements are taken to be Gaussian 
distributed random numbers with zero mean and variances 

(hln) = (Cn) = ^ n. (11) 

The random tunnelling interaction is taken to have zero mean and variance 



( 12 ) 


We assume that Ej lies in the middle of the N ND states, that is, Ej = 0. Following Refs, 
we introduce a spreading width Tj through the golden rule 


12.13 


Tj — 27rnj. 


(13) 


We maintain doubts about the meaningfulness of Eq. 
as a width 


^ ) regarding its interpretation 

T7[]. For practical purposes, however, it is a change of variable from vj to 


r^- 


All quantities of dimensions energy are to be understood to have units of Dj, where Dj 
denotes the mean spacing of the ND states around Ej. Thus the En (and the Eq and Vdq in 
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Eq. ([TSD) which are generated from the DGOE are to be understood to have been unfolded 
such that the and Eq have mean spacing equal to unity. Thus we may write Wjj{E) as 


Wjj{E) 




27r E Ef 


N 

E 




+ ^ /2 


(14) 


where the gn, n = 1, are Gaussian distributed random numbers with zero mean and 
unit variance. 

Model B is represented by the matrix 


( Ej 

Vjd 

0 \ 

n 


0 

0 \ 

Vjd 

Ed 

Vdq 


0 

r7 

J- N 

0 

V 0 

^dq Eqdq'q J 


1 0 

0 rl6q>qj 


q = -1. 


(15) 


We assume here that \J) couples to only one state of normal deformation; \d), which has 
energy Ed, with strength Vjd- This special state is subsequently mixed with other ND 
conhgurations with energies Eg by a residual interaction Vdg- Now Wjj{E) becomes 


N 


Wjj{E) = [VjdY E 


[cd{n)\^ 


■^^E-En + i{Vl + I) /2 


(16) 


where Cd{n) denotes component d of the nth eigenvector of the sub-matrix of the hrst term 
of Eq. ([T5|) obtained by excluding the hrst row and the hrst column. Now the Eq are 
eigenvalues of a random Hamiltonian h which is real symmetric and whose matrix elements 
are taken to be Gaussian distributed random numbers with zero mean and variances 

{Kq) = 9- (17) 


The residual interaction is also taken to have zero mean and variance 



E 

N' 


We put Ed = Ej = ^. Introducing 




we can write 


Eq. (p^ 


as 


Wjj(E) 


Nlcjnf 


(18) 


(19) 


( 20 ) 


Thus comparing Eq. (p!4[) with Eq. (p0|) we see that, although the meaning of Tj is diherent 
for the two models, the diherence between model A and model B boils down how much the 


distribution [cd{n)] dihers from that of 
dramatic A dependence (see Fig. (1) in Ref. 


This diherence is not trivial as [cd{n)\ 
]I5|] ). The inclusion of the factor T in Eq 


has a 

■ m 

makes clear that model A and model B are only comparable when Vj^ is of the order iVuj. 
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Model A is precisely equivalent to that of Refs. when A = 1. The real part of Eq. 

]) used in model B is equivalent to what is used in Refs. however, we calculate 

PM 


the average intensity integrated over the energy: FY', whereas Refs, 
tunnelling probability which is more closely related to Wjj{E). 


calculate a 


Note that F^' can be written as 


FT 


TA /•'» 

= ^ / dE- 
2tT J-oc 


-^72 


(B - B, - A‘,(E))2 + 4-(l + 


( 21 ) 




where 


A'j(E) = ReM'jXB), 


( 22 ) 


and 


r^(E) = -2lmWjj{E). 


(23) 


Ignoring the shift Aj{E) altogether and assuming that the width Tj{E) has the energy 

J 




independent value Tj one obtains the principal result of Ref. [12| that 


FT 








(24) 


We now present numerical calculations of F^' with N = 50, T} = O.OlDj. An ensemble 
average was performed over 100 realisations in Fig. (Q) and over 1000 realisations in Fig (^. 
The effect of increasing T^, identical to that obtained by increasing the energy averaging 
interval I, is to broaden rj(F) (Eq. (^)) and thus push F^' closer to the approximation 
given by Eq. (^). This is in line with what it is reported in [M who obtain Eq. (E^ in 

r'f 

the limit > cx) for their two level model. 

In our calculations we put + I = 3Dj. With this choice one may describe what Ref. 
calls the overlapping resonance region. Ref. |jT^ gives the impression that the relative 


intensity is independent of T}. Whilst we agree with that the ratios bj = and 
are of principal importance in understanding the decay out, it can be seen from Eq. (0) 
that FY' is only independent of T} if Aj{E) and rj(E) are constant. 

Fig. I shows the dependence of FY' on A, the strength of the mixing amongst the ND 
states, for several values of bj, for both model A and model B. For model A the variation of 
FY’ with A is rather slight compared to model B. This is because the A dependence of model 
A is contained in the eigenvalues En which are unfolded to have unit mean spacing. Model 
B has a further and more signihcant A dependence contained in the eigenvectors Q(n). For 
model B, FY' decreases with decreasing A to a value which is limited by the value of + /. 
Note that FY' can change at most by a factor of about 5 by varying A. 

Fig. ^ shows FY' as function of bj for some values of A, calculated using model B. 
The calculations for model A are not shown as they can barely be distinguished from the 
calculation for A = 1 using model B. The effect of changing A is to move the value of bj (and 
hence J) at which the decay out occurs. Thus from Fig. we conclude that the decay out 
is slightly hindered by increasing A. 
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Regarding Refs. |1I^,p!5| which report an increase in the tunnelling probability of several 
orders of magnitude with increasing A we do not consider ourselves at odds with this work 
since, as already mentioned above, we do not calculate the same quantity. A further differ¬ 
ence between model B of this paper and Refs. is that their author places Id) at the 

position ^ thus making the distribution Cd{'n) asymmetric. This would correspond in our 
calculation to making the difference Ej — Ed non-zero (we see no reason not to set Ej = Ed). 

An investigation of the roles of and I we postpone to a subsequent paper. The 


indicate that Eq. (|^, is important when is a small fraction of 


results of Ref 
Dj (non-overlapping resonance region). 


that Ef'- falls 


It was already found from phenomenological analysis some years ago 
exponentially with decreasing spin. We conclude here that the chaotic nature of the ND 
states, as classihed by A, cannot account for such behaviour. The exponential drop in the 
intra-band intensity must be due to the spin dependence of the tunnelling matrix element 
contained in bj. The calculation of bj is not trivial and we refer the reader to [^,^, which 
continue the work of Refs. @-|TI[ , for some recent calculations. 

A.J.S. thanks J.A. Tostevin for his comments on an early version of this paper. This 
work was supported by FAPESP. 
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Fj (dimensionless) 



X (dimensionless) 


FIG. 2. Calculated attenuation factor Ff'' as a function of A for some values of bj 
thin lines where calculated using model A and the thick lines using model B. 













(dimensionless) 



p-L 

FIG. 3. Attenuation factor F^' as a function of 6j = p4 for some values of A, calculated using 
model B. 
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